Mean field solution of the Ising model on a Barabasi- Albert network. 
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The mean field solution of the Ising model on a Barabasi- Albert scale- free network with ferromag- 
netic coupling between linked spins is presented. The critical temperature T c for the ferromagnetic 
to paramagnetic phase transition ( Curie temperature) is infinite and the effective critical tempera- 
ture for a finite size system increases as the logarithm of the system size in agreement with recent 
numerical results of Aleksiejuk, Holyst and Stauffer. 
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The Ising model is the starting point for the study of 
second order phase transitions and cooperative phenom- 
ena. Traditionally it has been solved exactly on periodic 
lattices in one and two dimensions |Q or in a Bethe lattice 
0. Recently attention has been addressed on the com- 
plex structure of scale-free networks || with power-law 
connectivity distribution P(k) ~ fc -7 . These networks 
describe for example biological systems where proteins 
are nodes and physical protein interactions are the links 
H . Moreover they constitute an intriguing substrate for 
the spreading of infectious diseases 0), and for percola- 
tion phenomena M . The Ising model has been applied to 
the evolving scale-free network with 7 = 3 [the Barabasi- 
Albert (BA) network §] by Aleksiejuk et al. (AHS) g 
and to random graph with arbitrary degree distribution 
[plfloj. By introducing a spin on each node of the net- 
work, they found that ordering occurs for temperatures 
T below the effective critical temperature T c , and that 
the nature of the phase transition is very different from 
that on periodic lattices. 

In this work I show that a mean field solution of the 
Ising model on the BA network can be treated as a Mat- 
tis model Jj"T| with the substitution £j — > ki and that 
the scaling of the effective critical temperature T c for fer- 
romagnetic transition increases as the logarithm of the 
system size in agreement with recent numerical calcula- 
tions ||. 

Let us consider a BA network of N nodes. Start- 
ing from a small number of nodes no and links mo 
(rto,77io << AT), the network is constructed iteratively 
by the constant addition of nodes with m links. The new 
links are attached preferentially to well connected nodes 
in such a way that at time tj the probability pij that the 
new node j is linked to node i with connectivity fci(ij) is 
given by 



fcj (tj ) 
i — — • 



(1) 



If N is large we can approximate the total number of 
edges present in the network at time tj, given by the 



sum Y2a=i with 2mtj and substituting into ([!]) the 
dynamic solution for the connectivity M 




we obtain 



Pi 



(2) 



(3) 



The adjacency elements of the network gjj are equal one 
if there is a link between node i and j and zero other- 
wise. Consequently their mean over many copies of the 
BA network has a tensor structure: 

hj] = Pij 

771 1 



2 ^Jtitj 



(4) 



The Ising- type spins Sj = ±1 placed on the nodes of a 
BA network are subjected to the Hamiltonian 



H 



(5) 



1,3 



with the local magnetic field hi and the coupling Jjj non 
zero only for nodes i and j connected by a link: 



Ji.j — J ^i,j 



(6) 



were J > 0, and e^j is the adjacency matrix of the BA 
network. 

The exact solution of the Hamiltonian is given by 



< Si > = < 



tanh[^(^ JijSj + hi)] 



> 



(7) 



The mean field equation for the mean local magneti- 
zation < Si > is given by 



A' 



< Si >= tanh[/3(J^[e y ] < Sj > +hi)\ (8) 

3=1 



1 



where we approximated < tanh(a;) > with tanh(< x >) 
and we performed the mean of the adjacency matrix [eij] 
over the different realization of the network. If we define 
S by 



I N 

i—1 



(9) 



after substituting [ei.j] by (|j) we can rewrite Eq.(| 



< Si >= tanh\P(JkiS + hi)}. 



(10) 
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Multiplying both sides of the equation by k i /2mN and 
summing over i we can solve (nQ) self-consistently for S, 



1 N 

i=l 



(11) 



Thus the BA Ising model looks like a Mattis model ]Tl| ] 
with the substitution 



(12) 



This is interesting because the Mattis model is a simple 
solvable model of spin-glass |L2| and suggest that in a 
spin-glass on a scale-free network there will be a noise 
term deriving both from the geometry of the system and 
the random interactions. 

The susceptibility- The susceptibility of the order pa- 
rameter S 1 , describing the response to a uniform magnetic 
field h, is given by 



E^i(i- < si > 2 ) 



8S_ _ 

dh ~ 2mN" (1 - ^/3J£ 4 fc?(l- < s t > 2 )) ' 



If we define the mean magnetization M as 



(13) 



N ^ 



< s t >, 



(14) 



FIG. 1. The analytical prediction for the ferromagnetic ef- 
fective critical temperature (solid line) is compared with the 
AHS simulations (open circles) with m = 5 as a function of 
the network size. 

The effective critical temperature- Using S as the order 
parameter of the system, and its self-consistent definition 
Eq.([ll]), it is easy to find the effective critical tempera- 
ture (T c ). As S ~ we can approximate the tanh(.T) 
with x and (3 with (3 C giving, 



2mN 
m 



N 



dt'l3 c JSk 2 (t') 



-Jp c Sln(N), 



i.e. 



To 
J 



ln(iV). 



(16) 



(17) 



This shows that for finite size system there is an effec- 
tive critical temperature that increases linearly with the 
interaction J and logarithmically with the number N of 
the nodes of the network. This result is in agreement 
with the numerical results of AHS || as it is shown in 
Fig. 1. In Fig. 1 we have compared the analytical results 
Eq. ( p"7j ) for the case m = 5 with the numerical results of 
AHS . However we know that a finite system cannot have 
a phase transition, therefore the true critical temperature 
goes to infinite in the thermodynamic limit N — > oo. This 
shows that in a scale-free network the ordered phase is 
the only allowed phase in the thermodynamic limit. 
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the susceptibility \ — ^j- is given by 



1 81 

= ^ {1 - <S * >2)[PJh dh +(3] 
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